Quantum Dynamics of a Driven Correlated System Coupled to Phonons 
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Nonequilibrium interplay between charge, spin and lattice degrees of freedom on a square lattice 
is studied for a single charge carrier doped in the t-J-Holstein model. In the presence of a static 
electric field we calculate the quasistationary state. With increasing electron-phonon (e-ph) coupling 
the carrier mobility decreases, however, we find increased steady state current due to e-ph coupling 
in the regime of negative differential resistance. We explore the distribution of absorbed energy 
between the spin and the phonon subsystem. For model parameters as relevant for cuprates, the 
majority of the gained energy flows into the spin subsystem. 

PACS numbers: 71.27.+a, 72.10.Di, 87.15. ht 



Introduction. — One of the outstanding contemporary 
challenges in condensed matter physics is to understand 
dynamics of interacting quantum systems exposed to 
an external perturbation. Advanced pump and probe 
techniques with few femtosecond time-resolution and 
broad-band THz spectroscopy were developed [u]^ to 
drive the system out of equilibrium and measure its 
nonequilibrium physical properties. These measurements 
were complemented by time-resolved photoemission spec- 
troscopy 15^6|, electron cristallography [7| and Raman 
scattering 8[. In the systems with competing interac- 
tions the most demanding task is to disentangle different 
elementary excitations arising at comparable energy-time 
scales. Possibly the most studied example of such ma- 
terials are cuprate superconductors where the influence 
of the electron-phonon (e-ph) interaction [9| on ultrafast 
dynamics was recently investigat ed, and different mecha- 
nisms were proposed [Sl, |5|, |8|, llCHl3j . Considerable effort 
has been devoted as well to study response of pump- 
excited Mott insulators [ij] and to unravel complex ther- 
malization processes of strongly correlated systems [l5| . 

Despite a considerable ongoing effort to understand 
nonequilibrium dynamics of quantum systems, a siz- 
able gap perseveres between theory and experiments. A 
theoretical insight into many-body quantum phenomena 
far from equilibrium has been obtained, among others, 
through works on nonlinear transport of half-filled Hub- 
bard systems 16l4l9j and Falicov-Kimball model [20| . It 
was found that to obtain a non-zero steady current in 
the metallic Hubbard model driven by a static electric 
field, it is crucial to couple the system to a heat reser- 
voir 18|, ll9| , which prevents an uncontrollable heating of 



the system. In contrast, the aim of our paper is to study 
a doped strongly correlated system coupled to phonons, 
where the energy gained by the motion of a charge carrier 
along the field is absorbed by quantum spin and phonon 
degrees of freedom which are all explicitly included in 
the model, and the influence of environment is supposed 



to show up on a much longer timescales. We study a 
single charge carrier doped into two-dimensional (2D) 
plane within the i-J-Holstein model, which is a proto- 
type model for the description of competing interactions 
in cuprates. Hereby, we address a fundamental, yet unre- 
solved question concerning the interplay between strong 
correlations and e-ph interaction in a driven quantum 
system far from equilibrium. 

So far, nonequilibrium response of the generalized 
Hubbard-Holstein model has been ajialyzed on a ID 
chain 
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and on 8-site 2D cluster [4|. While qualita- 
tively different behavior is expected in ID systems due 
to a spin-charge separation |22| . a detailed investigation 
of 2D systems is still pending. By generalizing recently 
developed method 23|,|24| , we drive the system by a static 
electric field until it reaches a quasistationary (QS) state 
with a constant current and a constant energy fiow to 
the system. To our knowledge this is the first study of a 
2D strongly correlated system (SCS) coupled to phonons 
where QS conditions are achieved. 

In this Letter, we explore two important aspects of 
nonequilibrium carrier dynamics: (i) We establish the in- 
fluence of e-ph coupling on the nonlinear transport prop- 
erties of a carrier in SCS. We show that the coupling to 
phonons decreases carrier mobility, however, it leads to 
an enhancement of QS current in the regime of negative 
differential resistance; (ii) We compare the energy ab- 
sorbed by the spin subsystem and the one absorbed by 
lattice vibrations. Taking into account model parameters 
fltting cuprates we find that the spin subsystem absorbs 
the energy from the field more efficiently than the lattice. 

Model and numerical method — We define the time- 
dependent i- J-Holstein Hamiltonian as 






A,'^ + H-c. 









(1) 



where Cj,a = cj,ct(1 — '^j.-o-) is a projected fermion op- 
erator, (/)ij (i) is time-dependent magnetic flux and (Ij) 
denote nearest neighbors. The charge carrier is cou- 
pled to Einstein phonons with the energy ujq via e-ph 
coupling constant g, where aj is the phonon annihila- 
tion operator and rij = X^o-^j.'^- ^^^ ground state at 
kg = (7r/2,7r/2) is calculated by exact diagonalization 



J=0.3, ^=0.2, CO =0.5 
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defined over a limited functional space (EDLFS) 

27|. To construct functions of the Hilbert space we use 

the off-diagonal parts of Eq. ([T]) in the basis genera- 

tor {l^f")} = [i/to(0ij = 0) + Hj + ifg]^H</'o>, where 
\Vo) = Ck(,|Neel) represents a translationally invariant 
state of a carrier in the Neel background. We switch on 
the static uniform electric field F along the diagonal at 
time t = and perform the time evolution by iterative 
Lanczos method [28|. Accordingly, we define the charge 
current along the diagonal j{t) and set (f>\j{t) = ~Ftj\/2 
for (positive) i,- and y-direction. We measure F in units 
of [io/eoa] and set to = eo = a = 1- Recently, the EDLFS 
method was reported to effectively calculate the QS state 
of a doped charge carrier within the 2D t~J model 1231 
as well as of the Holstein polaron 



2J1. The strength 



of the numerical method is in construction of the Hilbert 
space that enables not only an accurate description of the 
ground state of the spin-lattice polaron, but it allows for 
enough extra spin and phonon excitations to absorb en- 
ergy, emitted by the field driven carrier, until the system 
reaches the QS state. 

Results. — We focus mostly on weak and moderate val- 
ues of e-ph coupling A = g'^ /StoUJo and different regimes 
of ujo, while keeping J = 0.3 constant. Results describ- 
ing the real-time propagation of a spin-lattice polaron 
are shown in Figllja)-(c) for A = 0.2 and wq = 0.5. After 
a short transient regime i/i_B ^ 1 the system enters the 
QS state with the steady current j{t) = j and the linear 
increase of the total energy, satisfying AE{t) = Fj (com- 
pare Fig.UJa) and (b)), where AE{t) = {H{t)) - {H{0)). 
The longest time of propagation in the QS regime is lim- 
ited due to a finite number of spin and phonon excita- 
tions in the Hilbert space acting as reservoirs for the en- 
ergy absorption. They are determined roughly by N^, 
and choosing N^ = 10 enables computation of QS quan- 
tities with high accuracy (see the inset of Fig. [Iljd) for 
comparison with lower N^). In this work, we are inter- 
ested in values of F when the response of the system is 
dissipative, i.e., j ^ 0. The time evolution for F <ti 1 
becomes adiabatic with j ~ 0, as pointed out in discus- 
sion of Ref. 23|. In the opposite limit F ^ 1, F should 
not exceed the threshold electric field Fth required for the 
breakdown of Mott insulator at half-filling. Estimating 
Fth from the recent DMFT study {Fth "^ 10 in our units, 
J = 0.3) [ITI, the values of F indeed fulfill this condition. 

We next turn to the calculation of energy flow to the 
spin and phonon subsystem. In contrast to the noncqui- 
librium studies of closed systems at half-filling where the 




FIG. 1: (Color online) (a) AE{t) and (b) j{t) for F = 
1.6, 2.8 and 4.0. We measure time in units of Bloch time 
ts, where cub ~ F/^/2. Dashed lines in (a) represent ex- 
trapolation of linear increase of AE{t) in the QS state and 
dashed lines in (b) the corresponding QS current j. We 
set J = 0.3, A = 0.2 and cjo = 0.5. Numerical accuracy 
of the time propagation is checked by the total-energy-gain 
sum rule AE{t) ^ F J^ j{t')dt' . (c) AEj{t) and AE^^oit) for 
F = 2.8. Dashed lines represent their extrapolation to the 
QS state, (d) Energy distribution ratio in the QS state f; 
vs. F. Inset: Convergence of AE{t) for Nh = 7, 8, 9, 10 with 
Nst = 3.7 X 10'*, 1.6x 10^6.7x 10^ and 2.9x 10^ respectively. 



current resp onse is considerably influenced by the Joule 
heating [29[ , the problem of a single carrier in dissipative 
medium enables investigation of the steady growth of en- 
ergy due to carrier propagation in initially undistorted 
background at T = 0. In Figdfc) we show AEj{t) and 
AEa}„ (t), i.e., expectation values of the second and fourth 
term of Eq. ([1]). The energy flows to both subsystems are 
determined by Pj.uio{t) = AEj^^g{t), defining the distri- 
bution ratio ij{t) — Puio{t)/'Pj{t)- In the QS state both, 
AEj{t) and Ai?^Q(i), reveal a linear time dependence 
(dashed lines in Fig.[IJc)) and therefore T]{t) = fj. 

In Fig.[l](d) we show fj{F) which exhibits only tiny vari- 
ation around the constant value. This result is rather 
surprising since F strongly influences the energy flow 
into the system, nevertheless, fj remains fairly field- 
independent. 

This result facilitates the investigation of the efficiency 
of the energy absorption through spin and phonon chan- 
nel when e-ph coupling and phonon frequencies are var- 
ied. With increasing ojq, fj increases at fixed A as seen 
in Fig. mja). We are particularly interested in the case 
fj = 1, i.e., when a propagating carrier deposits equal 
amount of the gained energy to the spin and the phonon 
subsystem alike. For this purpose we calculate wo,i)=i 
and plot it vs. A in Fig. [21(b). Remarkably, a scaling 
i^o.fj=i ^ 1/vA is found. In the weak coupling regime 
the energy flow to phonons dominates over the flow to the 
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FIG. 2: (Color online) (a) f\ vs. a;o for different A and J = 0.3. 
At Wo = 0.2, ?7 is calculated for A < 0.2, while for higher loq 
a A^h-independent f\ is obtained until A = 0.4. (b) Red dots 
denote values of ciJo,f;=i for which the condition f; = 1 is ful- 
filled, and dashed line represents a fit a;o,J7=i = c/v^A, with 
c = 1.83. Left and right side of the plot correspond to the 
strong coupling (SC) and weak coupling (WC) regime of e-ph 
coupling, respectively. WC-SC crossover was calculated by 
the EDLFS method pal and is in agreement with other nu- 
merical methods pJ0l.l3ll|. Filled ellipse represents parameters 
as relevant for cunrates. 
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FIG. 3: (Color online) j-F characteristics of the QS state, 
(a) j-F at loq — 0.5 and (b) at A = 0.2. Solid line shows j-F 
at A = 0. Inset of (a): ^ — j/F (squares) calculated from the 
dashed lines in (a), LR mobility jj,lr ~ 2ttD (triangles) with 
charge stiffness D as calculated in Ref . [33| . Maximal current 
shows slight dependence on A^^ (compare j-F at A = with 
Ref. [23]), and contributes to discrepancy between fi and (j-lr- 
Inset of (b): crossover field Fo between PDR and NDR regime 
vs. ujo for different A. Dashed line shows Fo for A = 0. 



spin subsystem only for large values of wq > 1- Recent 
photo-emission 26|, |30| and optical experiments 32l |33| 
on cupratcs that were interpreted within the i- J-Holstcin 
model, assigned the realistic A to be in the interval 
[0.2,0.3] and Wq ^ 0.2, as indicated by the filled ellipse 
in Fig. mb). In the parameter regime as relevant for 
cuprates, the majority of the absorbed energy via the 
charge carrier driven by the constant electric field flows 
in the spin subsystem. 

Wc next explore the influence of e-ph coupling on car- 



rier's QS propagation. In Fig. [3] we plot j-F character- 
istics calculated in the QS state. Solid line represents 
j-F characteristics for a plain t-J model at J = 0.3, as 
calculated in Ref. [23[- In the latter work, it was shown 
that a regime of positive differential resistivity (PDR) 
at small F evolves into a negative differential resistivity 
(NDR) regime at crossover field Fq ~ 2.3J. The effect 
of increasing A, displayed in Fig. |3l[a) at ojq = 0.5, is 
to extend the region of PDR toward larger F, and to 
decrease the carrier mobility ^ as shown in the inset. 
A similar tendency is observed when increasing wq while 
keeping A fixed, see Fig.|3l^b). The carrier mobility in the 
PDR regime calculated from the j-F characteristics is in 
qualitative agreement with the linear response (LR) the- 
ory [231 , and allows one to detect the limit of LR regime 
occurring below the crossover field _Fo- The variation of 
_Fo with A and ujq is shown in the inset of Fig. [3jb) . 

While a decrease of mobility with larger e-ph coupling 
in the PDR regime can be intuitively understood due to 
increased scattering on phonon excitations, the most in- 
triguing result of Fig. [3] represents the phonon-induced 
enhancement of j in the NDR regime. In general, ap- 
pearance of NDR regime for large F is a consequence of 
limited degrees of freedom contained in the model that 
are available to absorb the excess energy, which impedes 
the carrier motion along the field. For instance, a NDR 
regime in 2D t-J model is characterized by pronounced 
transverse oscillations of the carrier that serve to emit 
the excess energy (gained by hopping along the filed di- 
rection) to spin excitations [23|. Alternatively, carrier 
propagation can be for large F interpreted on the basis of 
Wannier-Stark (WS) states where phonon assisted hop- 
ping between these states leads to nonzero j [2J, [SJ] . In 
this picture one can explain two phenomena observed for 
large F in Fig.[3lja): (i) with increasing F^ the overlap be- 
tween neighboring WS states mediated by the combined 
phonon- and magnon- assisted hopping decreases, hence 
j decreases with F, (ii) increasing the e-ph interaction 
seems to boost the already existing magnon- mediated 
overlap between WS states at A = 0. In contrast to low 
F, where additional scattering on phonons at A > di- 
minishes the carrier mobility, at large F the increase of 
j is due to opening of additional channels for deposit- 
ing the excess energy through simultaneous emission of 
magnons as well as phonons. This explains the seemingly 
counterintuitive effect observed in our results, i.e., that 
the current is enhanced due to an increased e-ph interac- 
tion. In a similar fashion, increasing of ujq enhances the 
phonon assisted hopping between WS states at fixed A 
and leads to an enhancement of j shown in Fig. ^h). 

Discussion and Conclusion. — By investigating the 
nonequilibrium charge dynamics in strongly correlated 
medium coupled to phonons, we showed that the ratio of 
the energy fiow to the phonon relative to the spin subsys- 
tem remains field- independent. Therefore, even though 
the values of F considered here are rather large, we ex- 



pect that the calculated distribution of the energy flow 
holds generally for any values of F for a weakly doped 
system where the dominant interactions are described by 
Eq. [T] In particular, for model parameters relevant for 
cuprates, the energy flow into the spin subsystem remains 
dominant. This result may signal stronger coupling of 
charge carrier to the spin system in comparison to e-ph 
coupling. An intuitive physical picture emerges when 
considering hopping of the hole in the spin background, 
simultaneously coupled to phonons. The hole can hop 
a few lattice spacings without exciting a single phonon, 
however each hop of the hole through a Neel background 
unavoidably generates spin excitations. 

The influence of phonons on carrier's QS propagation 
shows strong dependence on the particular regime of the 
j-F characteristics. While increasing e-ph coupling leads 
to a decrease of the carrier mobility at small F, the QS 
current increases due to increased coupling to phonons 
for large F, where the system enters NDR. Even though 
the calculation of a stable QS current in recent theoreti- 
cal studies of SCS is pos sible only by coupling the system 
to the heat reservoir [18|, |l9|, |29| , a qualitative connection 
can be made with our study where, in similar manner, 
the role of quantum phonons is to absorb the excess en- 
ergy. For instance, when the coupling to the heat reser- 
voir of a driven Hubbard model [l9| is increased, a shift 
of the maximal current in the j-F characteristics is ob- 
served accompanied with the increase of j in the NDR 
regime. These features are consistent with the influence 
of phonons in Fig. [31 Therefore, this work represents an 
example how to stabilize a QS state of a driven system 
where the work done by the field is absorbed exclusively 
by quantum degrees of freedom within the model. 

Our investigations of carrier dynamics under the static 
electric field are limited to modeling driving-induced low- 
energy intraband excitations. Focusing on the time- 
independent driving does not require consideration of 
interband transitions as long as the field, used in our 
calculations, remains below the threshold value for the 
dielectric breakdown of a Mott insulator. To mimic the 
situation realized in ultrafast experiments where high- 
frequency pulses are used, one has to combine mecha- 
nisms emerging due to both intraband as well as inter- 
band transitions. Still, a purely quantum mechanical de- 
scription of noncquilibrium dynamics in correlated multi- 
band models coupled to phonons remains an outstanding 
theoretical challenge in the field of driven systems. 
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